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ABSTRACT 


Several  empirical  models  for  predicting*  the  metabolic  response  to  a  lowered 
body  temperature  have  been  evaluaced  against  available  data  of  young  healthy  males 
immersed  in  cold  water  under  resting  conditions.  Nude  immersions  took  place  in  20 
and  24°  C  water  for  1  h  and  clothed  immersions  took  place  in  10  and  15°  C  for  3 
h.  The  data  were  pooled  according  to  low  and  high  percent  body  fat  (%BF). 
Decreases  in  the  mean  weighted  skin  temperature  (Tg^)  ranged  from  5.3  to  11.0°  C 
and  decreases  in  the  core  temperature  (T  )  ranged  from  0.56  to  1.54°  C,  while 
increases  in  the  metabolic  rate  over  the  immersion  period  ranged  from  34  to  256 
W.  Through  regression  analysist  an  inverse  relationship  between  %BF  and  the 
metabolic  response  for  a  given  lowered  and  lowered  T£  was  established.  When 
this  relationship  was  explicitly  applied  to  the  models,  significant  improvements  in 
their  predictive  capability  were  found.  Variables  such  as  body  weight,  body  surface 
area  and  the  rate  of  change  of  were  not  found  to  contribute  to  the  predictive 
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INTRODUCTION 

Several  studies  (/.~F>)  have  demonstrated  that  ■  iividual  of  low  subcutaneous 
body  fat  (BZ‘)  cools  more  quickly  than  one  of  high  BF  when  immersed  in  cold 
water.  This  increased  rate  of  cooling  is  associated  with  an  increased  metabolic 
response  in  the  form  of  shivering.  Shivering  is  the  body’s  method  of  generating 
heat  to  reduce  or  halt  the  rate  of  cooling.  It  is  natural,  therefore,  to  assume  that 
an  increase  in  metabolic  rate  correlates  to  a  decrease  in  deep  body  temperature. 
Indeed,  this  has  been  the  basis  of  several  models  that  predict  the  metabolic 
response  for  exposure  to  cold. 

Thermoregulatory  signals  that  determine  the  metabolic  response  are  thought  to 
be  generated  by  thermoreceptors  located  in  the  skin  and  the  hypothalamus.  In 
practice,  a  mean  weighted  skin  temperature  (T^)  is  used  to  define  the  signal  from 
the  skin.  To  approximate  the  signal  from  the  hypothclamus,  the  temperature  of  a 
particular  core  site  such  as  the  esophagus  ox  rectum  is  used.  In  the  present  study, 
the  temperature  of  these  core  sites  will  be  represented  by  Tc. 

The  models  that  predict  metabolic  response  from  measurements  of  T.k  and  Tc 
range  from  a  simple  linear  dependence  upon  these  two  variables  (6,7)  to  a  more 
complex  dependence  based  upon  the  differences  between  measured  and  fixed  or  set- 
point  temperatures  (8-11).  Most  models,  however,  do  not  diffarentiate  the  response 
for  individuals  of  varying  BF.  That  is,  the  predicted  metabolic  response  is  often  a 


function  of  only 


and  Tc,  yet,  for  a  given  lowered  and  lowered  T£,  there 


is  evidence  that  an  individual  of  low  BF  has  a  higher  response  than  an  individual 


of  high  BP.  This  stems  from  the  early  cold  water  immersion  studies  of  Gee  and 
Goldman  (12)  and  Bynum  and  Goldman  (13)  which  were  later  analysed  by  Strong 
et  al.  (7)  and  from  the  data  of  McArdle  et  al.  (14)  and  Toner  et  al.  (15)  used  in 
the  present  study.  Strong  et  al.  (7)  recognized  the  different  metabolic  responses 


between  individuals  of  varying  morphology  and  proposed  separate  model  equations 
f w.  small-lean,  average  and  heavy-fat  individuals. 

The  present  study  does  not  address  how  the  various  thermoregulatory  signals 
are  generated  or  integrated  to  produce  the  observed  metabolic  responses';  nor  does  it 
address  whether,  and  if  so  how,  these  properties  vary  between  individuals  of 
different  morphology.  Instead,  this  study  examines  how  well  the  above  referenced 
models  predict  the  metabolic  response  for  a  given  lowered  and  lowered  Xi. 

,  Using  the  more  recent  data  of  McArdle  et  al.  (14)  and  Toner  et  al.  (15),  the 
predictive  capability  of  these  models  are  compared  through  regression  analysis.  It  is 
then  demonstrated  that  the  predictive  capability  of  these  models  can  be  Lmprovdd 
by  explicitly  accounting  for  BF  through  an  inverse  relationship  between  the 
metabolic  response  and  BF.  Finally,  a  residual  analysis  of  the  best  model 
predictions  is  presented  to  examine  whether  the  assumption  of  normality  has  ni>t 
been  violated  and  whether  the  dependencies  on  Tsk  and  Tc  have  been  accounted 
for  adequately. 

METHODS 

Data.  Data  were  obtained  from  both  nude  (n~8)  and  clothed  (n=10)  yousg 
healthy  males  immersed  in  cold  water  under  resting  conditions.  Nude  (subjects 
wore  a  nylon  swim  suit)  immersions  took  place  in  20  and  24°  C  water  for  1  h  (14) 
and  clothed  immersions  took  place  in  10  and  15°  C  water  for  3  h  (15).  Clothing 
consisted  of  a  neoprene  dry  suit  (3  mm  thickness)  and  one  of  two  polypropylene 
undergarments  with  an  average  total  insulative  value  (based  on  calibration  with  in 
immersed  manikin)  of  0.95  (0.15)  and  1.28  clo  (0.20  ma  °C/W),  respectively. 
Rectal  and  skin  temperatures  were  measured  with  a  rectal  probe  and  waterproofed 
thermal  sensors.  Metabolic  rates  (MR)  were  determined  by  standard  techniques  of 
open-circuit  spirometry. 
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All  data  were  separated  according  to  low  (<14%)  and  high  BF  (individual 

%BF  ranged  fro*r>  *3  to  23.1%)  and  grouped  according  to  exposure.  T’  '  amber 

of  subjects  per  group  varied  between  4  and  5.  LBF  and  EBF  are  used  to 

designate  the  low  and  high  BF  groups,  respectively.  Although  each  clothed  subject 

was  exposed  to  the  two  immersion  temperatures  of  10  and  15°  C,  not  all  completed 

the  targeted  duration  of  3  h  and,  in  some  cases,  data  were  incomplete  v**)-  In 

these  circumstances,  the  number  of  subjects  in  the  group  was  reduced  from  5  to  4. 

Group  averages  were  adjusted  accordingly  when  these  data  were  applied  to  the 

models.  For  the  LBF,  these  average  values  ranged  as:  %BF  from  9.3  to  10.0,  body 

weight  (WT)  from  69.8  to  72.4  kg,  body  surface  area  (A^)  from  1.86  to  1.93  m® , 

decreases  in  from  5.3  to  11.9°  C,  decreases  in  Tc  from  0.98  to  1.54°  C,  and 

mean  increases  in  the  metabolic  rate  above  the  pre-immersion  resting  value  (AMR) 

from  52  to  250  W.  For  the  HBF,  average  values  ranged  as:  %BF  from  16.8  to 

19.7,  WT  from  74.3  to  80.0  kg,  A^  from  1.85  to  1.97  ma ,  decreases  in  from 

6.4  to  11.4°  C,  decreases  in  T  from  0.56  to  1.09°  C,  and  AMR  from  34  to  86  W. 

c 

Ratee  of  change  with  time  of  Tg^  and  T£  were  appioximated  by  linear 
interpolation. 

Models.  The  models  examined  are  listed  in  Table  1  and  although  this  list  is  not 
exhaustive,  it  does  represent  a  wide  cross-section  of  available  models.  The  first 
three  models  predict  MR  and  use  the  rectal  temperature  to  represent  T£  while  the 
remaining  predict  AMR  and  use  the  head  core  temperature  to  represent  T£.  The 
original  parameter  values  (Pi,  P2,  etc.)  were  converted,  if  necessary,  to  conform  to 
metabolic  rates  in  W  and  are  list  d  in  Tables  2  and  S. 


The  STRONG  model  (7)  was  originally  applied  to  groups  of  low  and  high  BF 
with  separate  parameter  values  (see  Table  3).  It  represents  the  simplest  model 
possible  with  and  Tc>  and  is  closely  related  to  the  TIMBAL  model  (6).  The 
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original  form  of  the;  TIMBAL  model  was  re-arranged  to  eliminate  unnecessary 
parameters  (with  respect  to  recession  analysis).  It  is  the  only  model  listed  to  use 
the  subject’s  Ap  and  the  rate  of  change  of  mean  skin  temperature  (!“.).  The 
HAYWARD  model  (S)  is  the  only  one  listed  to  use  the  subject1^  WT.  The 
inclusion  of  and  WT  into  these  models  will  be  examined  for  their  affect  on  the 
predictive  capability  of  the  models. 

The  WISSLER  model  (11)  expanded  upon  the  HAYWARD  model  by  including 
the  subject’s  rate  of  change  of  core  temperature  ft ).  The  form  shown  in  Table  1 
is  a  simplified  version  of  the  original  model  and  is  applicable  when  AMR  changes 
slowly,  as  was  the  case  with  the  present  data.  The  original  model  also  included  a 
term  with  ;  however,  this  term  decays  exponentially  with  a  half-time  of  30  s 
when  l1 ,  >  -1.5°C/min,  which  the  data  we  used  did  not  include,  so  this  term  was 
intentionally  left  out. 

The  NADEL  model  (8)  is  similar  to  the  HAYWARD  model  except  for  an 
additional  term  involving  only  T.k-  This  term  was  added  (8)  to  describe  the 
observed  metabolic  response  without  the  necessity  of  changing  the  hypothalamic  set- 
point  temperature  or  P3  in  the  me  del.  The  STOLWUK  model  (10)  emphasized  the 
contribution  of  this  additional  term  by  introducing  an  exponent  of  2.  Both  models 
consider  the  parameters  P2  and  P3  as  physiological  set-points,  to  be  identified  with 
temperatures  of  "standard"  man  under  a  condition  of  thermal  neutrality. 

The  PRESENT  model  is  based  on  the  STOLWIJK  model  with  the 


modification  that  the  set-point  values  and  Tfi  were  assigned  the  average 

of  the  measured  pre-immersion  values  and  remained  fixed.  For  example,  T.k,  ,set  ot 
a  particular  group  was  the  average  of  all  pre-immersion  values  of  T.k  for  that 
group.  This  procedure  allowed  the  use  of  "personalized"  set-points  which  did  not 
necessarily  match  the  values  of  "standard”  man.  For  data  used  in  this  study, 
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Tg^  ge^  ranged  from  32.79  to  33.50°  C  and  Tc  g^  ranged  from  37.14  to  37.39°  C; 
these  ranges  fell  outside  the  values  3  '7  and  36.96°  C,  respectively,  used  for 

"standard"  man  (10). 

Since  the  data  we  used  provided  rectal  temperature  as  the  only  Pleasured  core 
temperature,  this  temperature  was  applied  in  place  of  the  head  core  temperature  of 
the  WISSLER,  NADEL,  STOLWIJK  and  PRESENT  models.  In  these 
circumstances,  the  functional  form  of  these  models  and  not  necessarily  their  origin** 
parameter  values  were  tested.  In  anticipation  that  an  inverse  relationship  between 
BF  and  the  prediction  of  the  metabolic  response  for  a  given  lowered  T.k  and 
lowered  Tc  exists,  all  models  were  modified  "a  priori"  by  dividing  by  (%BF)X  where 
x  is  an  estimated  parameter.  These  modified  versions  are  designated  by  the  affix 
BF-modified. 

Analysis.  A  regression  analysis  was  performed  on  all  models  using  the  BMDP 
statistical  package  (16).  Significant  (a  =  0.05)  improvement  in  the  prediction  of  t 
particular  model  when  parameter  values  changed  or  when  a  parameter  was  added  to 
the  model  was  tested  using  Fisher’s  F-test  as  outlined  by  Mekjavic  and  Morrison 
(17>.  Significance  between  predictions  of  the  models  listed  in  Table  1  could  not  be 
tested  since  the.  models  were  fundamentally  different  from  one  another.  However,  in 
these  cases,  predictions  were  compared  by  the  models'  sum  of  squared  residuals 
(SSR)  given  by 


n=l 


where  N  is  the  number  of  cases  (grouped  observations),  is  the  case  weight 
(equal  to  1  for  groups  of  4  subjects  and  1.25  for  groups  of  5  subjects)  and  eQ  is 
the  residual  or  difference  between  the  observed  and  predicted  value.  In  essence,  the 
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SSR  is  a  measure  of  how  well  the  model  prediction  fits  the  data.  Models  with  the 
lowest  SSR  values  were  deemed  better  predictors  th  v  those  of  higher  SSR  values, 
although  no  significance  to  these  comparisons  could  be  given. 

The  residuals  of  the  model  that  predicted  best  overall  were  analyzed 
graphically  as  outlined  by  Draper  and  Smith  (It').  For  convenience,  the  "unit 
normal  deviate”  form  (e^/s)  of  the  residuals  was  plotted  where  s  is  the  square  root 
of  the  residual  mean  square  given  by 

s  =  ySSR/df  (0) 

and  df  is  the  number  of  degrees  of  freedom  (equal  to  N  -  p  where  p  is  the  number 
of  parameters  estimated). 

To  establish  whether  a  BF  dependency  in  the  prediction  of  MR  and  AMR 
exists,  each  model  in  their  unmodified  form  was  applied  to  the  LBF  and  HBF 
groups  separately,  similar  to  the  way  the  STRONG  model  was  originally  applied. 
Because  the  variability  of  BF  within  these  subgroups  is  small,  any  dependency  of 
the  metabolic  response  upon  BF,  if  it  exists,  is  effectively  removed.  Therefore,  the 
resultant  SSR  should  then  approximate  the  lowest  possible  value  without  introducing 
any  new  parameters.  To  compare  these  SSR  values  with  those  obtained  by  treating 
the  combined  data  (LBF+HBF)  when  using  the  same  unmodified  model,  the  sum  of 
both  groups,  i.e.,  SSR(LBF)+SSR(HBF),  was  used  so  that  all  the  residuals  were 
accounted  for.  If  little  difference  between  these  two  values  was  found,  then  it 


would  be  concluded  that  BF  is  not  a  critical  variable  in  the  predation  of  metabolic 
response  for  a  given  lowered  and  lowered  Tc-  However,  if  a  large  difference 
was  found,  then  BF  would  be  considered  a  critical  variable  and  its  role  could  then 


oe  quantified  through  modifications  of  the  models  such  as  described  by  the  BF- 
modified  versions.  An  appropriate  version  is  one  in  which  using  the  combined  data 
results  in  a  value  of  SSR(LBF-f HBF)  close  to  the  sum  of  SSR(LBF)-f SSR(HBF) 
using  the  separated  data  as  described  above. 
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RESULTS 

Comparison  Between  Original  and  Regressed  Estimates  for  ComL'ied  Data.  Table  2 
lists  the  model’s  (a)  original  (where  applicable)  and  (b)  repressed  parameter 
estimates  and  their  corresponding  SSR  when  applied  to  the  combined  data  of  the 
LBF  and  HBF  groups.  The  predictions  improved  significantly  for  a  ]  models  (where 
original  parameter  estimates  were  given)  using  the  regressed  parameter  estimates. 
The  following  comments  pertain  to  the  regressed  versions.  The  subscripted  alphabet 
in  parenthesis  affixed  to  SSR  refers  to  the  corresponding  set  of  parameter  estimates 
listed  in  Tables  2  and  3.  Among  the  regressed  versions,  the  WISSLER  and 
PRESENT  models  predicted  best  and  the  HAYWARD  model  predicted  worse  (see 
SSR^j  values  plotted  in  Fig.l). 

To  check  the  role  of  Aq  in  the  prediction  of  MR  in  the  TIMBAL  model,  this 
variable  was  removed  and  the  subsequent  regression  led  to  a  SSR  value  of  127,773. 
This  value  is  smaller  than  the  value  133,416  obtained  by  leaving  A^  in  the  model. 
Thus,  using  the  present  data,  the  inclusion  of  AD  worsened  the  predictive  capability 
of  the  TIMBAL  model.  A  similar  result  was  found  when  the  variable  WT  was 
removed  from  the  HAYWARD  model.  In  this  case,  the  resultant  SSR  was  123,146, 
considerably  smaller  than  the  value  158,933  obtained  by  leaving  WT  in  the  model. 

As  a  further  cheek,  the  variables  A^  and  WT  were  introduced  independently 
as  factors  to  the  PRESENT  model  equation  to  examine  whether  the  inclusion  of 
these  variables  could  improve  this  model's  predictive  capability.  The  subsequent 
regression  by  including  AD  led  to  a  SSR  value  of  34,694  and  that  obtained  by 
including  WT  was  05,112.  Both  values  exceeded  the  SSR  value  of  83,150  obtained 
when  neither  of  these  variables  was  used,  thus  no  improvement  was  achieved  by 
their  inclusion  into  the  model. 
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A  very  high  correlation  (r  >  0.00)  was  found  between  Pi  and  P4  of  the 
TIMBAL  model,  suggesting  that  *,k  w“  an  unnecessary  variable  to  »r  "  \de  for 
the  present  data.  This  is  further  substantiated  by  noting  that  the  SSR  value 
obtained  by  removing  Ajj  from  the  model  was  not  significantly  different  from  the 
v*lve  of  128,0x0  obtained  using  the  STRONG  model  which  is  identical  in  form 
exape  for  the  inclusion  ;*  the  term. 

A  linear  depeniv' ...  .vns  found  between  the  parameters  of  the  NADEL  model. 

,  This  was  not  unexpected  since  the  model  equation  (see  Eq.o)  could  have  been 
rewritten  without  the  second  term  such  that  P3(new)=P3+P4/Pl.  This  simplified 
form  was  not  used  by  Nadel  et  a!.  (8)  for  reasons  stated  earlier  but  was  proposed 
by  Montgomery  (19),  and  when  applied  to  the  regression  analysis,  it  was  found  to 
yield  the  same  SSR.  A  marginal  improvement  over  the  NADEL  model  was 
obtained  using  the  STOLWUK  model.  Because  of  the  exponent  of  2  raised  on  the 
second  term  of  this  model,  no  linear  dependency  of  the  parameters  was  found. 
Comparison  Between  Using  Combined  Data  and  Separated  Data.  Table  3  lists  the 
model’s  (a)  original  (applicable  only  in  the  case  of  the  STRONG  model)  and  (d) 
regressed  parameter  estimates  and  their  corresponding  SSR  when  applied  to  the  data 
of  the  LBF  and  HBF  groups  separately.  The  prediction  of  th*  STRONG  model 
was  improved  significantly  using  the  regressed  parameter  estimates  over  the  original 
parameter  values.  The  SSR  values  of  the  HBF  group  were  consistently  lower  than 
those  for  the  LBF  group  for  all  models  indicating  a  greater  variability  in  the 
metabolic  response  for  the  LBF  groups. 

To  compare  these  results  with  those  obtained  by  treating  the  combined  data, 
the  sums  of  SSR^(LBF)  +  SSR^(HBF)  are  plotted  alongside  SSR^j(LBF-f  HBF) 
in  Fig.l.  In  all  cases,  the  former  sum  is  subtaatiaily  smaller.  This  result  confirm* 
that  the  predicted  metabolic  response  to  a  lowered  T.k  and  lowered  has  a  BF 
dependency. 
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Functional  Role  of  BF.  Dividing  each  model  by  (%BF)X  (i.e.,  BF-modified  version) 
aignif  ’y  improved  the  prediction  over  the  unmodified  version*  using  the 

combined  data  in  all  cases.  Furthermore,  the  SSR(LBF-fHBF)  values  obtained  this 
way  were  close  to  the  sum  of  SSR^j(LBF)  +  SSR^(HBF)  indiciting  that  an 
inverse  relationship  of  BF  to  metabolic  response  was  appropriate  to  consider.  The 
regressed  estimates  of  x  varied  among  the  models,  but  for  practical  purposes, 
rounded  values  are  presented.  These  are  0.5  for  the  models  that  pred  .  MR  and 
1.0  for  the  models  that  predict  AMR.  The  regressed  parameter  estirr  .tes  of  the 
BF-modified  models  using  these  exponents  and  the  corresponding  SST  ^BF-rHBF) 
are  listed  in  Table  2.  These  SSR  values  are  also  shown  in  i.  Rounding-off 
the  exponent  x  was  possible  because  of  en  asymptotic  standard  deviation  of  about 
10%  in  the  estimate  of  x  and  because  of  a  high  correlation  between  x  and  some  of 
the  other  model  parameters. 

Residual  Analysis  of  the  PRESENT  BP-Modified  Model.  A  residual  analysis  of  the 
PRESENT  BF-modified  model  was  conducted  since  this  model  provided  the  lowest 
SSR.  Figure  2  shows  the  normalised  residuals  plotted  against  the  a)  predicted 
values  of  AMR,  b)  time,  c)  and  d)  Tc<  Since  05%  of  the  normalised  residuals 
fell  within  the  limits  (-1.06,  1.06),  this  distribution  was  assumed  normal  (18).  In 
addition,  no  apparent  trend  is  evident  in  the  plot  of  residuals  against  the  predicted 
values  of  AMR;  therefore,  the  residuals  were  also  considered  randomly  scattered. 


When  the  normalised  residuals  of  the  nude  and  clothed  groups  plotted  against 
time  in  Fig.  2(b)  are  considered  separately  (recall  that  the  nude  immersions  lasted  1 
h  and  the  clothing  immersions  lasted  3  h),  the  residuals  are  randomly  scattered; 
that  is,  no  time  dependency  is  apparent.  The  plot  of  residuals  against  shown 
in  Fig.  2(c)  reveeied  a  slight  tendency  from  overprediction  at  low  T.k  to 
underproduction  at  high  This  tendency  can,  if  desired,  be  reduced  without 
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comprcmising  the  randomness  of  the  other  plots  by  using  a  smaller  exponent  in 
-  place  of  2  for  the  <  nd  term  in  the  BF-modified  version  of  Eq.7.  The  final  plot, 
Fig.2(d),  of  residuals  against  T  indicated  random  scatter.  These  results  suggest 

V 

I 

that  anomalies  are  few  and  small;  hence  further  modification  of  the.  model  offers 
little  benefit  and  was  not  pursued. 

DISCUSSION 

The  results  suggest  that  an  individual’s  metabolic  response  to  a  lowered  T.k 
and  lowered  T  is  attenuated  by  his  %BF.  The  physiological  mechanism  underlying 
this  attenuation  is  beyond  the  scope  of  the  present  study.  However,  this 
attenuation  is  consistent  with  the  fact  that  increased  BF  promotes  thermal 
insulation  (1-3)  and,  therefore,  at  a  given  lowered  T»k  and  lowered  Tc,  less  heat 
production  is  required  to  maintain  deep  body  temperature. 

The  wide  disparity  between  the  original  parameter  values  and  the  regressed 
parameter  estimates  for  the  unmodified  version  of  all  models  may  be  explained,  in 
part,  by  differences  in  the  experimental  conditions  from  which  data  were  obtained  to 
derive  these  values.  Strong  et  al.  (7)  used  data  of  whole  body  immersion  in  water 
temperatures  between  20  and  36°  C.  Timbal  et  al.  (6)  and  Hayward  et  al.  (0)  used 
data  of  immersions  to  neck-level  in  water  temperatures  oi  15  and  26°  C,  and  10°  C, 
respectively,  and  Nadel  et  al.  (8)  used  data  from  exposure  to  cold  air.  Differences 
in  data  may  also  explain  why  the  regressed  parameter  estimates  proposed  by 
Mekjavic  and  Morrison  (17)  for  the  HAYWARD  "'odel  differed  widely  from  those 
we  found.  The  data  used  by  Mekjavic  and  Morrison  (17)  involved  immersions  of 
lean  individuals  in  10°  C  water.  These  authors  found  the  residuals,  when  plotted 
against  T»k'  displayed  a  shape  that  was  characteristic  of  the  average  discharge 
frequency  of  cold  receptors  in  the  skin.  This  bell-shape  with  a  maximum  at  about 
25°  C  is  not  apparent  in  Fig.2(c)  where  T#^  ranges  from  21  to  31°  C.  To  confirm 
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that  this  remit  was  model-independent,  a  residual  analysis  on  the  predictions  using 
the  original  HAYWARD  modT  \t  with  the  regressed  parameter  estimates  (set  (b) 

in  Table  2)  did  not  revsal  any  bell-shape  distribution  when  plotted  against  V 

«  ! 

The  addition  of  the  term  in  the  TIMBAL  model  did  not  improve  this 

model’s  predictive  capability  when  compared  to  the  results  of  the  STRONG  model 
which  is  similar  but  excludes  the  term.  We  ore  uncertain  about  the  value  of 
l1  used  in  the  WISSLER  model  despite  it’s  good  predictive  performance  because  of 
the  unusual  regressed  parameter  estimates.  It  is  interesting  that  the  TIMBAL  and 
HAYWARD  models  predicted  better  without  the  inclusion  of  the  Ap  and  WT 
variables,  respectively,  suggesting  that  these  variables  are  superfluous  in  the 
prediction  of  MR  using  the  present  data.  This  finding  was  further  supported  when 
Ap  and  WT  were  introduced  independently  in  the  PRESENT  model  and  its 
predictive  capability  worsened  in  both  cases. 

The  functional  form  of  the  HAYWARD  model  is  the  same  as  the  NADEL 
mo  '  c  asidering  that  the  second  term  of  the  NADEL  model  can  be  left  out 
without  affecting  its  predictive  capability,  yet  there  is  a  large  difference  between 
thek  SSR  values.  This  difference  is  attributed  partly  to  the  inclusion  of  WT  in 
the  HAYWARD  model  recalling  that  the  SSR  decreased  from  158,953  to  123,146 
when  WT  was  removed  from  the  model.  The  remaining  difference  is  largely 
attributed  to  the  prediction  of  MR  on  the  one  hand  and  AMR  on  the  other.  This 
is  because  the  variability  in  the  initial  value  of  MR  is  present  when  predicting  MR, 
as  in  the  case  of  the  HAYWARD  model,  but  this  variability  is  removed  when 
predicting  AMS  (recall  that  AMR  is  the  increase  in  metabolic  rate  above  the 
measured  pre- immersion  value).  Greater  accuracy  can  therefore  be  expected  from 
models  that  avoid  this  initial  variability.  This  is  further  supported  by  the  relatively 
poor  predictions  of  the  other  two  models,  STRONG  and  TIMBAL  that  also  predict 
MR. 


12 


Thu  difference  in  how  the  metabolic  response  is  predicted  also  helps  to  explain 
why  the  square  root  inverse  relationship  7  works  best  for  models  predicting  MR 

and  why  the  unit  inverse  relationship  of  BF  works  for  models  predicting  AMR. 

1 

For  the  same  increase  in  metabolic  rate,  the  relative  change  in  AMR  la  larger  than 
for  MR.  It  follows,  therefore,  that  relative  differences  in  AMR  between  the  LBF 
and  HBF  groups  are  also  larger  than  relative  differences  in  MR  between  these 
subgroups.  These  differences  can  be  adequately  balanced  by  dividing  the  original 
models  that  predict  AMR  by  %BF  and  those  that  predict  MR  by  \[% BF. 
CONCLUSIONS 


Regressed  parameter  estimates  for  all  unmodified  models  improved  the 
prediction  of  metabolic  responses  significantly  over  tbe  original  parameuer  values. 
This  suggests  that  when  these  empirical  models  are  used,  the  values  of  the 
parameters  should  correspond  to  exposure  conditions  from  which  they  were  derived. 
In  addition,  since  the  esophageal  temperature  represents  more  closely  the  head  core 
temperature  than  does  the  rectal  temperature,  the  predictive  capabilities  of  the 
WISSLER,  NADEL,  STOLWUK  and  PRESENT  models  may  improve  with  the  use 
of  the  esophageal  temperature,  if  available. 

Applying  each  model  to  the  LBF  and  HBF  groups  separately  demonstrated 
that  BF  is  a  critical  variable  to  consider  when  predicting  the  metabolic  response  for 
a  given  lowered  T.k  and  lowered  Tc>  Applying  an  inverse  relationship  of  %BF  to 
the  metabolic  response  (i.e.,  the  BF-modified  versions)  significantly  improved  all  the 
models’  predictive  capability.  The  PRESENT  model  predicted  best  overall,  both  in 


the  unmodified  form  when  applied  to  the  LBF  and  HBF  groups  separately  and  in 
the  BF-modified  form  when  both  groups  were  combined.  The  best  alternative  if 
specific  thresholds  or  "personal"  set-points  were  unavailable  was  the  STOLWIJK 
model  using  the  regressed  parameter  estimates  given  in  Tables  2  and  3. 


13 


ACKNOWLEDGEMENTS 

The  authors  gratefully  acknowledge  Igor  B.  M  vie  and  Michael  M.  Toner 
for  valuable  discussions  during  the  course  of  this  work,  and  Edna  E.  Safran  for 
preparing  the  manuscript.  1 

The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
authors  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  official  documentation. 

Approved  for  public  release;  distribution  unlimited. 


14 


REFERENCES 

1.  r  sr  PB,  Daniels  F.  Relationship  between  skinfold  ness  and  body 
cooling  for  two  hours  at  15°  C.  J  Appl  Physiol  1956;8:409-416. 

2.  Carlson  LB,  Hsieh  CL,  Fullington  F,  Eisner  RW.  Immersion  hi  cold  water 
and  body  tissue  insulation.  Aviation  Med  1958;20:145-152. 

3.  Nadel  ER,  Holmer  I,  Bergh  U,  Astrand  P-O,  Stolwijk  JAJ.  Energy  exchanges 
of  swimming  man.  J  Appl  Physiol  1974;36(4):465U71. 

4.  Keatinge  WR.  The  effects  of  subcutaneous  fat  and  of  previous  exposure  to 
cold  on  the  body  temperature,  peripheral  blood  flow  and  metabolic  rate  in  cold 
water.  J  Physiol  1960;153:166-173. 

5.  Sloan  REG,  Keatinge  WR.  Cooling  rates  of  young  people  swimming  in  cold 
water.  J  Appl  Physiol  1973;35(3):371-375. 

0.  Timbal  J,  Loncle  M,  Boutelier  C.  Mathematical  model  of  man's  tolerance  to 
cold  using  morphological  factors.  Aviat  Space  Environ  Med  1976;4?(9):058-964. 

7.  Strong  LH,  Gee  GK,  Goldman  RF.  Metabolic  and  vasomotor  insulative 
responses  occurring  on  immersion  in  cold  water.  J  Appl  Physiol  1985;58(3):964-977. 

8.  Nadel  ER,  Horvath  SM,  Dawson  CA,  Tucker  A.  Sensitivity  to  central  and 
peripheral  thermal  stimulation  in  man.  J  Appl  Physiol  1970;29(5):6O3-6O9. 

9.  Hayward  J§,  Eckeraon  JD,  Colli*  ML.  Thermoregulatory  heat  production  in 
man:  prediction  equation  based  on  skin  and  core  temperatures.  J  Appl  Physiol 
1977 ;42(3)  :377-384. 

10.  Stolwijk  JAJ,  Hardy  JD.  Control  of  body  temperature.  In:  Lee  DHK,  ed. 
Handbook  of  Physiology.  Reactions  to  Environmental  Agents.  Am  Physiol  Soc. 
Betheada,  1077:45-68. 

11.  Wissler  EH.  Mathematical  simulation  of  human  thermal  behavior  using  whole 
body  models.  In:  Heat  transfer  in  medicine  and  biology  v.l.  Plenum  Press.  New 
York,  1985:325-373. 

12.  Gee  GK,  Goldman  RF.  Heat  loss  of  man  in  total  water  immersion. 
Physiologist  1973;16:318. 

18.  Bynum  GD,  Goldman  RF.  Whole  body  cooling  with  protective  clothing  during 
cold  water  immersion.  Physiologist  1974;17:191. 

14.  McArdle  WD,  Magel  JR,  Gergley  TJ,  Spina  RJ,  Toner  MM.  Thermal 
adjustment  to  cold-water  exposure  in  resting  men  and  women.  J  Appl  Physiol 
1084;56(6):  1565-1571. 

15.  Toner  MM,  Holden  WL,  Foley  ME,  Bogart  JE,  Pandolf  KB.  Influence  of 
clothing  and  body-fat  insulations  on  thermal  adjustments  to  cold-water  stress. 
Aviat.  Space  Environ.  Med.  (submitted). 


15 


16.  Dixon  WJ.  BMDP  statistical  software.  Univ  California  Press.  Berkeley, 
1083. 

17.  Mekjavic  IB,  Morrison  JB.  Evaluation  of  predictive  formulae  for  de (,«.’* mining 
metabolic  rate  during  cold  water  immersion.  Aviat  Space  Environ  Med  1986;57: 
671-680. 

t 

18.  Draper  NR,  Smith  H.  Applied  regression  analysis.  John  Wiley  and  Sons, 
New  York,  1966. 

19.  Montgomery  LD.  A  model  of  heat  transfer  in  immersed  man.  Annals 
Biomed  Eng  1974;2.l0-46. 


16 


T 


FIGURE  CAPTIONS 

I 

1.  Sum  of  squared  residuals  for  the  combined  data  using  the  regressed  parameter 
estimates  (SSR^(LBF+HBF)),  for  the  separated  data  using  the  regressed  parameter 
estimates  (SSR^(HBF)),  and  for  the  combined  data  using  the  BF-modified 
regressed  parameter  estimates  (SSR^  (LBF+HBF)). 


2.  Normalized  (unit  normal  deviate  form)  residuals  of  the  PRESENT  BF-modified 
model  plotted  against  (a)  predicted  increase  in  metabolic  rate,  (b)  time,  (c)  mean 
weighted  skin  temperature  and  (d)  core  temperature.  The  closed  circles  represent 
the  nude  condition  and  the  open  circles  represent  the  eletbed  condition. 


STRONG  TIMBAL  HAYWARD  WISSLER  NADEL  STOLWIJK  PRESENT 

MODEL 


MEAN  WEIGHTED  SKIM  TEMPERATURE  (°C)  CORE  TEMPERATURE 


i 


H 

*S 

i 

o 


■S5 

c 

re 

Je 

«-r 


i! 

<u 

5 

o 


4) 

V) 

il 

8 

fi£ 

o 

1 

a 

3 

•ij 

o 


O 

o' 


Jt 

M 


U1 


wt 


4) 

Wt 

Jt 

W) 


-K 

M 


a. 

■ 


n 

Q. 


n 

a. 

I 


o  , _ „ 

►?  m 

w  a. 


in 

a. 


a. 


£ 


i  t 

£  2 
w  + 


a 

+ 


CL 

+ 

'~C 

H 

I 


CD 

fc- 

« 

1 

w 

P2) 

o 

H, 

(P3 

(P3 

• 

£ 

m 

a. 

*- 

1 

« 

• 

Iff 

■ 

m 

a) 

-t3 

je 

CM 

H 

CN 

Q 

W» 

JC 

M 

1 

5 

W 

K 

CL 

• 

l 

#- 

i 

I 

« 

M 

m 

• 

T*1 

▼H 

J£ 

i/i 

CN 

CN 

ji 

vi 

u 

CM 

CL 

CL 

t- 

Ou 

CL 

H 

CL 

i 

*E 

c 

• 

• 

• 

• 

« 

• 

E 

UJ 

*-i 

CL 

Q 

< 

WT 

Pi 

pi 

pi 

pi 

II 

U 

II 

u 

u 

ii 

u 

C 

X3 

O' 

O' 

CX 

cc 

CX 

cx 

(X 

TO 

2 

2 

5 

H* 

<3 

<3 

< 

<3 

Z 

o 

'M 

< 

CQ 

cx 

< 

£ 

cx 

ID 

_l 

to 

ID 

ex 

>• 

to 

O 

h- 

< 

< 

</> 

1— 

X 

$ 

z 

* 


o 

h 

tn 


H 

Z 

ID 

(n 

ID 

CX 

CL 


T 


[• 


1 


.15 

c 

A 


Hu 


L 

* 

Hi 


WO  K 

"»C  « 

MIH 

K  K  © 

®  «  a 

•  a 

eel 

•-*  Y 

hhH 

W  C  WO 

so« 

•  ®  r» 

C  Y 

<aI 

o  « 

n  y  k 

«  3  •* 

CYC 

BH» 

•  as 

•x  a 

*0  tt 

DMA 

t#N 

WNH 

«  •  N 

a  tern 

t 

WK 

H 

KMC 
M  H 

<-MAC 

W»4 

9ttt 

C4 

•  y 

Y  M 

U5  ©  Y 
H  H 

CM 

t 

£1 


K  N  Y 


•  II 


•h  III 


I  I  I 


I  I 


«9« 


I  Cl 
I 


HtO 


I  I  I 


NtOl 
»  N 
I  I 


*  C4  < 
•  • 

•  •i 


•  i 


U. 

O 

* 


n\ 

C  Y 

c  k  w 

Y  M  C 

U»  10  • 

O  M  Y 

K 

* 

<M 

»  • 

•  •  • 

•  • 

*  •  • 

■  •  • 

•  •  • 

JO 

U 

10  *0 

H69 

-<#• 

t*  K  K 

c  c  c 

©  K  K 

t  1 

• 

WO  W 

<4  M  U> 

Y  M  V 

Y  M  M 

M  M  M 

M  M  M 

Y 

H 

to 

W* 

• 

« 

%- 

* 

Y 

s 

e 

a 

•r*. 

Me, 

K  M  «H 

Y 

• 

csteJ  (Mi 

w4  K 

®  £5  SJ 

•*  <■  © 

•  tN 

NO^ 

Acer 

M  Y 

t 

•.» 

•  • 

•  •  * 

«  •  « 

»  *  • 

•  •  * 

•  •  • 

•  « 

« 

•* 

O  K 

W»  Y  ^ 

\S 

<S4  54  r* 

«  Y  M 

X 

M 

M 

Y  Y  V 

Y  M'S* 

6$*$  M 

«  M  « 

4» 

• 

• 

A 

u, 

■ 

Y  K  K 

n 

*- 

•*  Y  -M 

MS 

Y  0W 

c  •  « 

M  Y  M 

fM  H 

C 

•  *•  V 

•  ••* 

KWH 

CMC 

•  • 

C4  • 

*-»| 

Y  H 

H41 

•  •  • 

•  •  « 

•  •  • 

«  •  « 

*  • 

1 

wj 

c*  » 

•  •• 

WON 

»»KH 

CM 

U» 

Y  © 

•H  K 

K 

*4 

X 

Y 

t 

• 

• 

A 

•— 

/N/N 

y-v^-v 

^\rs 

%- 

<Y  O 

*  -a  u 

•  .tf  U 

•  J3  U 

m  x  u 

•  JQ  U 

JB  0 

»— 

v*w 

'-~r^ 

ww 

Y 

© 

fi 


6 

o 

•s 

►- 

i 

a 

-  _j 

UJ 

M 

mm 

X 

•4 

< 

_» 

c 

UJ 

© 

1* 

M 

u 

_l 

• 

fc 

» 

h 

M 

Q 

O 

u 

*- 

H 

< 

H 

< 

H 

* 

© 

t- 

X 

© 

Z 

c 

■< 


\ 


*>  ‘  r  •» 


(Oltl 

^  N 

no 

0£ 

«  «  ©  •» 

©  rl 

N  'C 

W 

«MI« 

*  N 

ID  CM 

m 

«  «  «  » 
NOtHK 

iet»n 
C4  10 

«  * 

•H  M 
10  «H 

*  •> 

N  M 
tH 

ft  • 

H  H 

OW 

K  •* 

© 

e  a 

*  © 

ft  o 

K  © 

H  01 

f-  * 

ID  oT 

a  ©' 

MO* 

J  aT 

CM 

CM  r4 

M 

C4 

•te 

cb/) 

di 


i  i  i  i 


i  i 


cm  f- 
cm  © 

CM  fH 
I 


I  I 


I  I 


© 

10  ft 

'Cl 

11)1 

*  ft 

h  r* 

K  CM 

«  * 

0.1 

•  • 

•  • 

•  . 

•o 

m 

«  l£ 

i  i 

•  w 

CM  K 

«*• 

1  1 

ft  T4 

CM 

• 

| 

1  1 

£s 

XJ 

o 

0 

• 

4> 

o 

• 

t» 

© 

• 

e 

a 

CM  ■* 

• 

• 

**l 

•  <0  <■  N 

K  H 

ft 

ft  CM 

«  » 

CM  N 

(/> 

0.1 

♦  ♦  »  « 

•  t 

«  • 

•  • 

•  • 

HU)«H 

ID  *S 

H  H 

t*  r* 

40  ID 

K  N 

1  1 

m 

L 

h  ft  ft 

**  H 

t  t 

m  m 

ft  ft 

ft  ft 

a 


•  WHt  ft 


•* 

c»  © 

ID  ID 

01  ft 

f- 

o* « 

•  ••r 

ft  ID 

•  • 

ft 

H  H 

ft  f 

CM  ft 

ft  ft 

ft  ft 

• 

• 

Of 

*  K 

»  © 

W  H 

«  H 

to  «-* 

Of  M 

ID 

*  * 

ft  K 

•  m 

K  ID 

ID  © 

**1 

#-  » 

©  * 

t  K 

*  • 

•  • 

♦  ♦ 

•  • 

♦  • 

0.1 

H  *4 

•f 

m  m 

©  O 

B  •• 

•  •• 

© 

«4 

*4 

i 

i 


33 

33 

33 

^0  © 

33 

Ik  Ik 

fa.  U- 

Ik  (k 

iku. 

tklk 

Ik  Ik 

U  VL 

©  © 

©  © 

©  © 

©c 

c© 

m  © 

©  e; 

-1  X 

JS 

JX 

JZ 

JX 

JX 

JX 

X 

- 

p 

© 

H 

O 

.J 

•£ 

to) 

X 

X 

< 

< 

-J 

-4 

m 

bl 

o 

m 

BB 

tt 

w 

_j 

9C 

m 

>- 

tt 

Q 

o 

u 

»- 

w 

< 

M 

«* 

H 

© 

to 

t- 

X 

* 

X 

M 

4. 

